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Comment on "Degenerate Wannier The- 
ory for Multiple Ionization" 

Recent Letter of Pattard and Rost ||l|] suggests new 
threshold laws for the processes of break up of a par- 
ticle into several charged fragments, which, regretfully, 
is incorrect. Previous studies, starting from the classical 
paper by Wannier [^ derived fragmentation cross section 
(T in a power-law form 


a(e -> 0) 


(1) 


where e is the energy excess above the break up thresh- 
old. The primary task of the theory is to evaluate the 
threshold index /i for different multiparticle systems, a 
comprehensive bibliography could be found in Ref. [Q. 
Pattard and Rost |^,Q] suggest the threshold law of novel 
functional form, namely, with an extra logarithmic e- 
dependence: 


a(e ^ 0) 


|ln< 


(2) 


According to the originally idea of Wannier g , for the 
Coulomb interaction between the fragments the thresh- 
old law is defined by the motion along the 'potential ridge 
as the constituent particles fly apart. This motion is un- 
stable; survival on the ridge corresponds to the system 
fragmentation, whereas sliding from it means that the 
break up is not achieved. The instability of the motion 
is characterized by a set of Lyapunov exponents A^. The 
threshold index /i is expressed via suvn of all Lyapunov 
exponents PM]. In the present Comment we discuss the 
key result obtained by Pattard and Rost in ||I||, and re- 
peated in [Q, namely, the novel threshold law (||). 

The standard approach to derivation of threshold laws 
relies on the analysis of classical trajectories in the vicin- 
ity of the potential ridge where the equations of motion 
are linearized over coordinates transversal to the ridge. 
The system is described by a set of coupled linear equa- 
tion that are differential over the time variable t. The 
time-dependence of the coefficients is eliminated by intro- 
ducing the effective time r. The solutions of these equa- 
tions (/(r) depend on the effective time as exp(AiT) where 
the set of Xi eigenvalues is found by solving the charac- 
teristic equation. The threshold index ^ is expressed via 
Ai, see for details Ref. P,p[ and bibliography therein. 

Pattard and Rost |jl^,ff essentially retain the traditional 
framework, but argue that the situation is changed when 
the Lyapunov exponents happen to be degenerate. They 
refer to the standard mathematical textbook H saying 
that "If n eigenvalues are degenerate, the general solution 
contains additional terms T'"'exp(Ar), {k < n)". Further 
they relate the logarithm in the alleged coordinate time- 
dependence, 

q{T) = r'' exp(AT) ee exp(AT + fclnr) {k < n) , (3) 

to the logarithmic factor in the energy-dependence (0). 


This reasoning is invalid. Consider harmonic vibra- 
tions of a system with many degrees of freedom (for ex- 
ample, a polyatomic molecule) around an equilibrium po- 
sition. Some eigenfrequencies might be degenerate due to 
symmetry reasons, or accidentally. The argumentation 
by Pattard and Rost fully applies to this case; however, 
as universally known, the logarithmic solutions (0) in fact 
do not emerge [see Eq. (23.6) in Ref. Q and subsequent 
discussion]. The reason is clear: in the harmonic ap- 
proximation all the normal modes are fully decoupled and 
hence a character of the time-dependence in each mode 
does not depend on whether some other degenerate mode 
exists or not. 

This argument is directly related to the motion in the 
vicinity of the potential ridge. The problem is described 
by the Hamiltonian that is quadratic both in transversal 
coordinates and momenta, similar to the harmonic ap- 
proximation for a polyatomic molecule. The only differ- 
ence is that the equilibrium is unstable and some eigen- 
frequencies are complex- valued ||^. Obviously this fact 
does not influence the functional form of solutions of the 
same equations of motion which still reveal no logarithmic 
terms. The degenerate Lyapunov exponents are related 
to the different modes which are fully decoupled, there- 
fore the logarithmic solutions (^, being mathematically 
feasible for general linear differential equations, do not 
emerge in the physical applications concerned. 

As a summary, the threshold law is given by Eq. dl^) 
while the modification (0) is never valid. 
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